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Abstract 

We study a particular texture of the neutrino mass matrix generated in super- 
gravity with bilinear R-Parity violation. The relatively high value of tan (3 makes 
the one-loop contribution to the neutrino mass matrix as important as the tree-level 
one. The atmospheric angle is nearly maximal, and its deviation from maximal mix- 
ing is related to the smallness of the ratio between the solar and atmospheric mass 
scales. There is also a common origin for the small values of the solar and reactor 
angles, but the later is much smaller due the large mass ratio between the light- 
est two neutrinos. There is a high dependence of the neutrino mass differences on 
the scalar mass tuq and the gaugino mass Mi/ 2 , but a smaller one of the mixing 
angles on the same sugra parameters. Measurements of branching ratios for the 
neutralino decays can give important information on the parameters of the model. 
There are good prospects at a future Linear Collider for these measurements, but 
a more detailed analysis is necessary for the LHC. 
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1 Introduction 



With a number of experimental results in atmospheric, solar, reactor, and accelerator 
neutrino physics, it has been established that neutrinos have mass and oscillate 0. This 
is a very important result in its own but, in addition, it is the first direct experimental 
indication that the Standard Model (SM) needs to be modified jif. 

In the SM neutrinos are massless. One popular mechanism for the generation of neu- 
trino masses is the see-saw mechanism, where a right handed neutrino field with a very 
large mass is added to the SM Q. The resulting neutrino mass is inversely proportional 
to this large mass. Another interesting and predictive mechanism is the radiative gener- 
ation of neutrino masses and mixing in a supersymmetric model Q that violates lepton 
number and R- Parity [^| with bilinear terms in the superpotential. Phenomenological con- 
sequences of R-Parity violating supersymmetry are very distinct from R-Parity conserving 
models p. 

Bilinear R-Parity breaking is an interesting mechanism for the generation of neutrino 
masses and mixing angles due to its simplicity and predictability 0, B|- It is a simple 
extension of the Minimal Supersymmetric Standard Model (MSSM) which includes no 
new fields and no new interactions. It differs from the MSSM in a handful of bilinear 
terms that violate lepton number and R-Parity, which cannot be eliminated with field 
redefinitions Neutrino masses and mixing angles are calculable and agree with ex- 
perimental measurements 0, Ull. Motivations for BRpV are for example, models with 
spontaneously broken R-Parity [12|, and a model with an anomalous horizontal U(l) 
symmetry |l3j], where BRpV appears without trilinear R-Parity violation. 

ResultsTom SuperKamiokande Q on atmospheric neutrinos gave strong evidence of 
the — v T oscillation of the atmospheric neutrinos with maximal or nearly maximal 
mixing, and gave strong evidence against the small mixing angle solution of the solar 
neutrino problem. Results from the Sudbury Neutrino Observatory (SNO) and the Kam- 
LAND experiment have confirmed the large mixing angle solution of the solar neutrino 
problem, showing that more than a half of the electron-neutrinos produced at the sun 
oscillate into other flavours before reaching Earth Q. Results from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) show temperature differences within the microwave 
background radiation, which combined with results from large scale structure give a bound 
on the sum of the neutrino masses jl4j . Finally, evidence for neutrinoless double beta de- 
cay, if confirmed, would show the Majorana nature of neutrinos and the non-conservation 
of lepton numbers 0] . 

There are several analysis of these experimental results ^(J. The 3a allowed regions 
for the neutrino parameters in 0] are 

1.4 x 1(T 3 < Am 2 2 <3.3xl(T 3 eV 2 

7.2 x 1(T 5 < Am 2 21 < 9.1 x KT 5 eV 2 

0.52 < tan 2 fl 23 < 2.1 (1) 

0.30 < tan 2 # 12 < 0.61 

tan 2 #13 < 0.049 
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which we show for reference. 

In this article, we re-analyze the possibility of having BRpV in a supergravity scenario, 
in which the scalar masses and the gaugino masses are universal at the GUT scale. The 
electroweak symmetry is broken radiatively but, contrary to the MSSM, sneutrinos acquire 
vacuum expectation values as well as the Higgs bosons. We give up the possibility that 
the 6j and Bi parameters (one for each lepton and analogous to the /i and B term in the 
MSSM respectively) are universal at the GUT scale, because otherwise there is no good 
solution for the neutrino physics compatible with experiments. 

We found solutions that have not been discussed previously in the literature. These 
solutions are characterized by a large value of tan/5 and, therefore, the importance of 
one-loop contributions to the neutrino mass matrix is enhanced. 



2 Neutrino Mass at Tree Level 

The superpotential of our BRpV model differs from the MSSM by three terms which 
violate R-Parity and lepton number, 



W = Wmssm + €tLiH v 



(2) 



where q have units of mass. We complement them with related terms in the soft la- 
grangian, 

£ = £mssm + Bi€iLiH u (3) 

where £?; also have units of mass. The presence of these terms induce vacuum expectation 
values Vi for the sneutrinos, which are calculated minimizing the scalar potential. 

At tree level neutrino masses are generated via a low energy see-saw type mechanism. 
Neutrinos mix with neutralinos, and the MSSM neutralino mass matrix is expanded to a 
7x7 mass matrix for the neutral fermions 



M N = 



m 



m 




T l 



(4) 



Here, M. x o is the usual 4x4 neutralino mass matrix, and m is 



m 



h'vi \gv x 



\g'v2 \gv 2 o e 2 
lgv 3 e 3 



(5) 



which mixes the neutrinos with the neutralinos. The matrix M^r can be diagonalized by 
blocks, and the effective 3x3 neutrino mass matrix turns out to be equal to 



M(°) = 



! T M ig 2 + M 2 g' 2 
-m ■ M o ■ m = r 

x 4det(M x o) 



A? A X A 2 A X A 3 
A X A 2 A 2 2 A 2 A 3 
AiA 3 A 2 A 3 A 2 3 



(6) 



3 



where we have defined the parameters Aj = /iVi + e^, which are proportional to the 
sneutrino vacuum expectation values in the basis where the e terms are removed from the 
superpotential. 

This mass matrix can be diagonalized with the following two rotations. 



V (o) 



t 1 


1° 





cos 



(o) 

23 



\ ( COS Bis 



sm y 2 3 



(0) 
23 



- sin 0. 

COS 023 



J 



V Bill eg 





1 





sin 0® \ 



cos 6® 



(7) 



where the reactor mixing angle in terms of the alignment vector A is 

Ai 



tanflg 



(A 



AS 



1 5 

2 



and the atmospheric angle is 



tan 



(o) 

23 



V 



(8) 



(9) 



As we will see later, despite the fact that tree level contribution to the heavy neutrino 
mass dominates over all loops, there are other contributions to the neutrino mass matrix 
that cannot be neglected. For this reason, the above tree level formulas will not be enough 
to explain the results. 



3 Supergravity and BRpV 

In Sugra-BRpV the independent parameters are 

m , Mi/2 , A , tan/3 , sign(/x) , e { , A; , (10) 

where mo is the universal scalar mass, Mi/ 2 is the universal gaugino mass, and A is 
the universal trilinear coupling, valid at the GUT scale. In addition, tan/5 is the ratio 
between the Higgs vacuum expectation values, and sign(/x) is the sign of the higgsino 
mass parameter, both valid at the EWSB scale. Finally, q are the supersymmetric BRpV 
parameters in the superpotential, and Aj are the parameters depending on the sneutrino 
vacuum expectation values. 

We use the code SUSPECT jl8| to run the two loops RGE from the unification scale 
down to the weak scale. The electroweak symmetry breaking is analogous to the MSSM, 
with the difference that there are three extra vacuum expectation values corresponding 
to the sneutrino vev's These vev's are small and constitute a small perturbation to 
the MSSM EWSB. 

Despite the fact that sneutrino vev's are dependent quantities since they are calculated 
from the minimization of the scalar potential, we remove from the group of independent 
parameters the B^s in favour of Aj = fiVi + eiVd as indicated in eq. (fTUj) . because they are 
more useful in describing the neutrino physics. 
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Our analysis will be centered around the SPS1 scenario in Sugra from the Snowmass 
2001 benchmark scenarios flti ]. which is defined by 

m = 100 GeV, M 1/2 = 250 GeV, A = -100 GeV, tan/3 = 10, // > (11) 

This scenario is typical of Sugra, with a neutralino LSP with a mass m x o = 99 GeV, and 
a light neutral Higgs boson with a mass just above the experimental limit rrih = 114 GeV. 

In this context we find several solutions for neutrino physics which satisfy the ex- 
perimental constraints on the atmospheric and solar mass squared differences, the three 
mixing angles, and the mass parameter associated with neutrino-less double beta decay 
For illustrative purposes we single out the following 

e 1 = -0.0004, e 2 = 0.052, e 3 = 0.051 , 

A : = 0.022, A 2 = 0.0003, A 3 = 0.039 , (12) 

This solution is characterized by 

Am 2 32 = 2.7 x 10~ 3 eV 2 , Am 2 21 = 8.1 x 10~ 5 eV 2 , m ee = 0.0036 eV 

tan 2 9 23 = 0.72 , tan 2 12 = 0.54 , tan 2 13 = 0.0058 (13) 

which are well inside the experimentally allowed window in eq. We note that the 
random solution in eq. ()12j) is compatible with e\ = A 2 = 0, i.e., the neutrino parameters 
in eq. (fT3*j) are hardly changed with this replacement. 



4 Texture of Neutrino Mass Matrix 

Among the solutions to neutrino physics that we have found in our model, there are a few 
textures |22j] for the effective neutrino mass matrix. Our study case in eq. ()12|) belongs 
to the most frequent one, which is 



Mf, ff 



m 



A A 
a a 
A a 1 



(14) 



with a ~ 0.5 — 0.8, A ~ 0.1 — 0.3, and m ~ 0.02 — 0.04 eV. To understand how this texture 
works we expand the neutrino masses and mixing angles in powers of A. Keeping terms 
up to first order, the three neutrino masses are 



m 
m 



IJ 2 



m 



Am + C(A 2 ) 

V5a 2 - 2a + 1) + C(A 2 ) 
\m{\ + a + V5a 2 - 2a + 1) + C(A 2 ) 



(15) 



and the rotation matrix that diagonalizes the neutrino mass matrix, denoted Upmns, is 

1 Xsgm/m u , 2 Xcgm/m U3 

UrMXs= A (1 a) c„ ' -sg ' fO(A 2 ) (16) 

-A/(l-a) Sg Cg 
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with 

—2a 

tan 26 = . 17 

1 — a 

In the same approximation, the atmospheric, solar, and reactor angles are given by 

tan 2# 23 = 2a/ (1 -a) + C(A 2 ) 
tan 6 12 = \sQin/m U2 + 0(\ 2 ) 

sin ^13 = \c e m/m V3 + £>(A 2 ) (18) 
while the atmospheric and solar mass differences are 



Am 2 32 = m 2 (l + a)V5a 2 - 2a + 1 + £>(A 



2\ 



Amjj = \m 2 



1 + 3a 2 - (1 + a)V5a 2 -2a + l + C(A 2 ) . (19) 



As an example, consider a = 1/2 and m = 0.04 eV. We find Am^ = 3\/5m 2 /4 « 
2.7 x 10" 3 eV 2 , and Am^ = (7 - 3\/5)m 2 /8 « 5.8 x 10~ 5 eV 2 , both in agreement with 
experiments. The third parameter which in this approximation does not depend on the 
small parameter A is the atmospheric angle, obtaining tan 2 823 ~ 0.4 from eq. ()17j) . This 
value is at the lower end of the allowed region, nevertheless, taking a = 0.6 we obtain 
tan 2 823 ~ 0.5, which is in better agreement with experiments. 

The fact that a is smaller than unity implies that the atmospheric mixing is not 
maximal. In the limit a — *■ 1, the atmospheric mixing approaches maximality, but the 
atmospheric mass Am 2 tm — > 4m 2 which is too large if m = 0.04 eV, and the solar 
mass Am 2 ol — > which is too small. Decreasing tan 2 8 a tm via decreasing a will decrease 
the atmospheric mass scale and increase the solar one, both towards acceptable values. 
Therefore, the value of a relates these three neutrino parameters, such that the non- 
maximal value for the atmospheric angle is connected to the smallness of the ratio between 
the solar mass scale and the atmospheric one. 

The previous considerations are modified by the non zero value of A. In the approx- 
imation we are working, the solar and reactor angles are proportional to the parameter 
A, thus they are small quantities themselves. Nevertheless, the presence of m U2 in the de- 
nominator of tan# 12 as opposed to m U3 in the denominator of tan #13, makes the reactor 
angle much smaller than the solar angle. In the case a = 1/2 and A = 0.2 we find for the 
solar angle tan 2 812 = 0.3, which is in the lower part of the allowed region and compatible 
with experiments. For the reactor angle we find tan 2 #13 = 0.017 which is well below the 
experimental upper bound. We stress that we use the complete numerical calculation in 
the rest of the article, rather than these approximated formulas. 



5 One Loop Contributions 

All particles in the MSSM contribute to the renormalization of the neutralino/neutrino 
mass matrix. One of the most important contributions comes from the bottom-sbottom 
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loops. In the gauge eigenstate basis this contribution is (21 

N c m b 



An 



167r 2 AS b C b n b^0 



~ — i e i A j + e A) + I — 1 



where we have defined 



9 /l 

71 U %ai ~ a2 



2 
a L — — y= -t w ai 



g'M 2i i gM lt i 
a\ = — — — , a 2 



«3 



J u , 2 



(g 2 M 1 + g f2 M 2 



2A - 2A 4A 
The main contributions to eq. ()20|) can be understood as coming from the graph 



(20) 



(21) 



(221 




a 3 A j - ej/n 



a 3 Ai - ei/n 



Here, neutrinos (in the gauge eigenstate basis) mix with Higgsinos who in turn interact 
with the pair bottom-sbottom with a strength proportional to the corresponding Yukawa 
coupling. Full circles indicate the projection of the sbottom mass eigenstate into right and 
left sbottom, which contribute with a sin 9~ b and cos 6~ b respectively. Open circles indicate 
the projection of the neutrino field onto the higgsino, proportional to the small parameter 
G^Aj — ej//i. The quark propagator contributes with a factor m^, and summing over color 
gives the factor N c . Finally, we have in eq. (J2(Jj) 

AB~f b * = B (0,ml,m 2 b ) - 5 (0, m? 2 , r^) (23) 

The one-loop corrected neutrino mass matrix, in first approximation, has the general 
form 

Ally = AA { Aj + B(eiAj + SjA^ + Ce^j (24) 

since all loop contributions can be expanded in this way. The terms of higher order in A 
and e have been neglected. 

Considering the solutions to neutrino physics whose effective neutrino mass matrix has 
a texture of the form in eq. (|14|). and including contributions from all one- loop graphs, 
we extract the numerical value of the above parameters and find A pa 7eV/GeV 4 , B pa 
-0.5 eV/GeV 3 , and C pa 9 eV/GeV 2 . 

Of the three parameters only A gets a contribution at tree level, and we estimate 

A<n = * M \ + f M > MeW/Gey* (25) 
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Clearly, the tree level contribution to A dominates over all one-loop graphs. This is not 
true for B and C because this two parameters are entirely generated at one-loop. 

The contribution to A, B, and C from the bottom-sbottom loops can be read from 
eq. (J2UJ). In the squark sector we have mg = 492, mi = 538 GeV, and sin 26^ = 0.88, 
which implies 

<f = -^ s i c i h l AB ~o' b2 « 9 - 8 eV / GeV2 ( 26 ) 

This result is very close to the actual numerical value, and underlines the fact that the 
bottom-sbottom loops are very important in this particular scenario. 

Considering that the value for B, in the supergravity model we are working with, is 
much smaller than A and C, we might in first approximation neglect it in eq. (J24j) . In this 
case, for the neutrino solution in eq. (|12[) we obtain the following approximated neutrino 
mass matrix, 

AA\ AAiA 3 

-2 

This form is precisely the texture observed in eq. (jl4j) obtained from the numerical results. 
Therefore, the zero in the neutrino mass matrix is there because A2, e\ ~ and because 
B is very small compared with A and C. The three matrix elements of order A in eq. (|14|) 
are explained by the fact that Ai has a numerical value smaller than the other three 
relevant parameters, as can be seen from eq. (|12|) . Finally, the parameter a in eq. (|14|) 
is smaller than unity because AA\ and Ce\ are comparable and of the same sign and 
because e 2 ss e 3 . 



Mf/ / 



Ce\ Ce 2 e 3 
AA t A 3 Ce 2 e 3 AAJ + Ce 2 



(27) 



6 Numerical Results 

In this section we study numerical results on the neutrino mass matrix, neutrino mass 
differences and mixing angles. We center our studies in the supergravity benchmark 
given in eq. 1)110. although we also explore the behavior of the neutrino parameters in 
the m — M1/2 plane. We look for solutions to neutrino physics with different values of 
the BRpV parameters and Aj, but concentrate our attention in the particular solution 
given in eq. (fT2~j) . 

First, we consider the supergravity benchmark in eq. (|TT|) and randomly vary the 
BRpV parameters and Aj. We look for solutions satisfying experimental restrictions 
on neutrino parameters according to the 3a intervals in eq. and also according to a 
relaxation of those cuts given by: 

1.2 x 1(T 3 < Am 2 32 < 4.8 x l(T 3 eV 2 

0.43 < tan 2 9 23 < 2.3 (28) 
5.1 x 10~ 5 < Am 2 21 <19xlO~ 5 eV 2 

motivated by previous allowed regions and shown in order to compare the effect of the 
improved analysis of the experimental data. 
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Figure 1: Solutions to neutrino physics passing all experimental cuts described in the text, 
working within a particular supergravity benchmark. 



9 



-0.10 -0.05 0.00 0.05 
e 2 (GeV) 



0.10 




-0.10 -0.05 0.00 0.05 
e 3 (GeV) 



0.10 



co 5 
i 

o 

- 1 4t- 
X 



CM 
> 



3 — 



2 — 



i i i i | i i i i | i n i i 

(c); 

+ 
+ 

V,**"*"" m ^ m / 
I < I I I ' ' I I ' I I ' I ■ ' ■ ' 



-0.10 -0.05 0.00 0.05 0.10 
r2> 



Aj (GeV ) 



CO 

I 

O 

t-H 

02 
> 



a 
< 



5 
4h 
3 — 
2 — 
1 — 



i i r+i i i i i i n i i i 14-1 i 



/(<*) : 



...J*... 







-0.10 -0.05 0.00 0.05 

t3> 



0.10 



A 3 (GeV") 



Figure 2: Atmospheric mass squared difference as a function of the four relevant BRpV 
parameters for the reference scenario: €2, 63, Ai, and A3. 



Solutions satisfying the relaxed cuts given in eq. (|2*5jl are displayed as green crosses in 
Fig. [TJ over the plane formed by the absolute value of the vector e and the squared root 
of the absolute value of the alignment vector A, both quantities measured in GeV. Two 
distinctive regions are observed, with low and large values of |A|, with the low value of |A| 
solutions harder to obtain. When the stringent cuts are implemented we find solutions 
only in the region of large |A|, and we represent them as red squares. 

Since the tree-level neutrino mass matrix depends on Aj only, and one-loop corrections 
depends on both Aj and e iy although dominated by e iy the position of the solutions in the 
plane |ej v/s |A| is an indication of how important loop contributions are. We stress 
the fact, nevertheless, that increasing values of tan/3 (which we keep constant in this 
study) increase the importance of one-loop corrections, as observed in eq. (|2TIj) due to the 
presence of the Yukawa couplings. In our case, as we will confirm in the following figures, 
the one- loop contributions to the neutrino mass matrix are very important. 

In Fig. El we have the atmospheric mass squared difference as a function of the four 
BRpV parameters €2, €3, Ai, and A3. The neutrino mass matrix has the texture shown in 



10 



eq. ()14|) , which implies an atmospheric squared mass difference Am^ given approximately 
by eq. (|19|) . The parameters m and a in eq. (|14|) are m = A A 2 + Ce| and a = Ce 2 /m, as 
can be read from eq. (|27|) . The scale m is quadratic in the parameters A 3 and €3, since the 
dependence of A and C on A's and e's is weak. The dependence of the atmospheric mass 
is obtained by replacing these expressions in eq. (fTUj) . but when a ~ 1/2 the atmospheric 
scale can be approximated even further obtaining, 

Am 2 2 ss fV^A 2 + Ce^Ce 2 (29) 

explaining the quadratic dependency of Am^ on €2, 63 and A 3 , in frames (j2k), (J2b), and 
(j2t0 respectively, and the mild dependency on Ai (hidden in the neglected terms of order 
A 2 ), as can be observed in frame (j2fc). The dependence on Ai become strong at high 
values of this parameter because in that case neglected terms are no longer small. 

We note that using the tree level formulae in chapter 2, the atmospheric mass scale 
would be given by Am^ P3 (A^|A| 2 ) 2 f=3 0.3 x 10~ 3 eV 2 , highlighting the inadequacy of 
the tree level formula. On the contrary, the approximated expression in eq. (|29|) gives 
a value 2.8 x 1CT 3 eV 2 , which is much closer to the value in eq. (|T3)l found using the 
complete calculation. 

In Fig. El we plot the tangent squared of the atmospheric angle, tan 2 2 3 . Using eq. (fTHj). 
or directly from the mass matrix in eq. (|27p. we find that the atmospheric angle satisfy 

tM ' 2 ^M + c(|- t !) < 3 °) 

This relation implies that if €2 approaches zero, the atmospheric angle 623 — > 0. This 
behavior is confirmed in frame (jSk)- On the other hand, if 63 approaches zero then 
023 — * 7r /2 because Ce 2 , is larger than AA 2 , and this explains the divergence of tan 623 in 
frame (jUb)- 

In frame (JHfc) we see again the mild dependency of the atmospheric parameters on A2, 
in this case the atmospheric angle. If this parameter becomes very large though, neglected 
terms of the order A 2 become important. Finally, the dependency of the atmospheric angle 
on A3 in frame (01) can be understood also from eq. ()3U|) since clearly if 1 7V3 1 grows then 
tan 6*23 decreases. 

From eq. 0, the tree level atmospheric angle satisfy 



no (p) 2A 2 A 3 
tan 20^ = A 2_ A 2 ( 31 ) 



and this relation clearly misses all the influence of the one-loop graphs to neutrino mass 
matrix seen in eq. ()30j) . Numerically, the approximated formula in eq. ()3Uj) gives tan 2 023 ~ 
0.42, which is close to the value in eq. (|13|). On the contrary, the tree level formula implies 
tan 2 9$ w 0. 

In Fig. E] we plot the solar mass squared difference as a function of the BRpV parame- 
ters e 2 , e 3 , Ai, and A 3 . In the case of Am^ the neglected terms of order A 2 in eq. (fTU|) are 
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Figure 3: Atmospheric angle as a function of the four relevant BRpV parameters for the 
reference scenario: e 2 , e 3 , A l7 and A 3 . 
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Figure 4: Solar mass squared difference as a function of the four relevant BRp V parame- 
ters for the reference scenario: e 2 , e 3; A 1; and A 3 . 
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numerically more important than in the atmospheric case, therefore, predictions based on 
this approximation are less accurate. 

In frame (|3k) we see the dependence of the solar mass on e^- This behavior can 
be understood considering that the parameter a is proportional to e|, and when this 
parameter goes to zero, the solar mass difference approaches zero like a 2 , as seen from 
eq. ©. 

In frame (@Jd) we see how the solar mass difference depends on €3. If 63 — > then 
the eigenvalue Ce 2 , decouples and becomes the heaviest neutrino. Of the other two, one 
neutrino is massless, and the solar mass difference becomes equal to the second neutrino 
mass squared. A growing e 3 will mix the massless neutrino with the heaviest, increasing 
the lightest neutrino mass, therefore, decreasing the solar mass difference, as observed in 
frame 

The dependency of the solar mass on Ai and A 3 can be understood only if we go 
beyond the simple approximation in eq. (fEJj) . Terms of order A 2 introduce a dependency 
on Ai and A3 such that A approaches to zero when these last parameters go to zero, thus 
explaining the behavior shown in frame (j3t) and (jUi). 

In Fig. we have the tangent squared of the solar angle, tan 2 # 12 , as a function of the 
BRpV parameters €2, €3, A 1; and A 3 . Working in the texture given in eq. (|14|) . the solar 
angle according to eq. (fTHj) is approximately given by 

AA 2 

tan# 12 = sin 6 23 (32) 

m U2 

The dependency on Ai is explicit and comes from the small parameter A in eq. (|T4"|) . As 
we know from eqs. (fTHj) and (JTHJ), the dependency of 6*23 and m U2 on A x is weak. The 
behavior of the solar angle on Ai seen in frame (03 ) is thus understood. 

The solar angle as a function of e 2 can also be easily understood noting that the 
parameter a is directly proportional to e|. According to eqs. (fTB]) and (fTSj) the second 
neutrino mass approaches zero when a — > 0, explaining the divergence shown in frame 
(JHk). Note that sin #23 also approaches zero when e 2 — * 0, but slower. 

The divergence of tan # 12 when A 3 — > is harder to understand from the approximated 
expression in eq. (|32|) . so we go back to the effective neutrino mass matrix in eq. (|27|). If 
A3 approaches zero then the upper-left element of the matrix decouples with a mass AA 2 . 
On the other hand, the lower-right 2x2 sub-matrix has a zero eigenvalue, implying that 
9 12 — > 7r/2, and therefore explaining the divergence shown in frame (JUi). 

In Fig. El we have chosen the neutrino solution given by the BRpV parameters in 
eq. (fT2|) . and vary the scalar mass m and the gaugino mass Mi/ 2 , looking for solutions 
that satisfy all experimental cuts. In this case, sugra points satisfying the experimental 
restrictions on the neutrino parameters lie in the shaded region. Solutions are concentrated 
in a narrow band defined by Mi/ 2 ~ 230 — 260 GeV and m ~ — 400 GeV. We note that 
in BRpV the LSP need not to be the lightest neutralino, since it is not stable anyway. 
For this reason, the region close to m ~ is not ruled out. 

Smaller values of Mi/ 2 are not possible because the atmospheric and solar mass dif- 
ferences become too large. The allowed strip is, thus, limited from below by the curve 
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Figure 5: Tangent squared of the solar angle as a function of the four relevant BRpV 
parameters for the reference scenario: €2, 63, A 1; and A3. 
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mSUGRA: tanj8=10, A=-100 GeV, /^>0 




m GeV 

Figure 6: Region of parameter space in the plane mo — M1/2 where solutions to neutrino 
physics passing all the implemented experimental cuts are located. Contours of constant 
atmospheric mass difference and angle, and solar mass difference are displayed. 

Am^ = 9.1 x 10~ 5 eV 2 . The dependency on M1/2 is felt stronger by the tree level con- 
tribution to the parameter A, given in eq. (J25j) . There we see that A decreases when the 
gaugino mass M1/2 increases, implying that the atmospheric mass decreases with M1/2, 
as seen in eq. (|2"9"|) . In addition, the solar mass difference is proportional to the parame- 
ter m 2 , which in turn is proportional to A, thus, the solar mass also decreases with the 
gaugino mass. 

Higher values of the scalar mass mo are not allowed because the atmospheric angle 
becomes too small. The allowed strip is, therefore, limited from the right by the contour 
tan 2 823 = 0.52. We can understand this behavior in the following way: the parameter 
C decreases with increasing mo due to the Veltman's functions, and this in turn makes 
tan 2 623 to decrease with the scalar mass. High values of the scalar mass are also limited 
from above because the atmospheric mass becomes too large. This can be explained from 
eq. (|29|) considering that the parameter C decreases with increasing m . 

Higher values of M1/2 are not possible because the solar mass becomes too small, 
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Figure 7: Solutions to neutrino physics in the plane formed with the atmospheric and the 
solar mass differences. For the three different values of mo =100, 150, and 200 GeV, we 
vary the gaugino mass Mi/ 2 . 



therefore, the allowed stripe is limited from above by the line Am^ = 7.2 x 1CT 5 eV 2 . 
As we already mentioned, the solar mass difference is proportional to the parameter m 2 , 
which in turn is proportional to A, and we already know that A decreases with increasing 
gaugino mass My 2 - 

In Fig.[7|we see from another point of view the dependence of the solar and atmospheric 
mass differences on the scalar mass m , and the gaugino mass Mu2- In the plane formed 
by the atmospheric and solar mass differences we plot four curves defined by a constant 
value of the scalar mass m = 100, 200, 300, and 400 GeV, and vary the gaugino mass 
in its allowed region, which is indicated in the figure. We keep fixed the values of the 
BRpV parameters 6j and A, given eq. (|12|). The two neutrino mass differences are clearly 
proportional to each other highlighting their common origin represented by eq. ([24]) . where 
the parameter A is controlled by tree-level physics and the parameter C is controlled by 
one-loop physics, and where both are equally important. Fig.[7|can be understood further 
when seen in relation with Fig. El 
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Figure 8: Partial decay width of a neutralino into a W and a lepton, measured in units 
of distance. 



7 Collider Physics 

In our model, lepton number and R-Parity are not conserved. One important consequence 
is that the lightest supersymmetric particle (LSP) is not stable, and will decay into SM 
particles. Since it is not stable, the LSP needs not to be the lightest neutralino, and 
whatever it is, its decays can be used to prove the BRpV parameters and the neutrino 
properties [l0|. In the supergravity benchmark point considered here, the LSP is the 
lightest neutralino, with a mass m x o = 99 GeV. 

One of the interesting decay modes of the neutralino is Xi ~^ W^F, where I = e, /i, r. 
This decay is possible because the neutralino mixes with neutrinos which in turn couple 
to the pair Wl, and also because the charged leptons mix with charginos and they in 
turn couple to the pair x\W . For this reason, the relevant couplings in this decay are in 
general very dependent on and A*. 

In Fig. |H1 we plot the inverse of the partial decay width (multiplied by the velocity of 
light to convert it into a distance) as a function of the most relevant BRpV parameters. In 
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5^ ^X° f 4 ■ X° 9'' ^X° 

Figure 9: Feynman diagrams relevant for the production of two neutralinos at the LHC . 

frame (jBK) we see the inverse of T(x\ — > We) as a function of Ai. In fact, for all practical 
purposes, the decay rate into electrons depends only on Ai. Since in first approximation, 
the coupling is proportional to A 1; the inverse of the decay rate behaves like A]" 2 , and this 
is seen in the figure. The values of Ai are limited by the solar parameters. The inverse of 
the partial decay rate Xi ~^ We is of the order of 30 — 40 /im, and it's an important part 
of the total decay rate. 

In frame (JBb) we have the inverse of V(xi — > Wfi) as a function of A 2 , and similarly 
to the previous case, the decay rate into muons depends practically only on A 2 . In our 
reference model in eq. fjl2jl we have A 2 ~ 0, but values indicated in the figure are also 
compatible with neutrino physics. The coupling of the neutralino to W and muon is 
proportional to A 2 , so the inverse of the decay rate goes like A^ 2 , and that is observed in 
frame (JBJd). Depending on the value of A 2 , the partial decay length vary from millimeters 
to more than a hundred meters in the figure. Therefore, this partial decay rate contribute 
little to the total decay rate of the neutralino. 

The inverse of r(x? -> Wt) is plotted in frames (JEfc) and 01) as a function of Aq 
and €3 respectively. The dependence on A 3 is stronger and similarly to the previous cases 
it goes like A3 2 . The dependence on €3 is weaker, and the inverse decay rate increases 
with this parameter. The inverse decay rate is of the order of 8/im, making it the most 
important contribution to the total decay rate. Neglecting any other decay mode, the 
total decay rate is near 6/iin. The ratios of branching ratios for our benchmark point in 
eq. (fT2|) are given by 

B{X\ ±VW _ 5 9 x 1(r5 gfrg ^We) _ 

B(x ^Wr)- 5 - 9Xl ° ' B(x°^Wr)-°- 32 (33) 

We note that if we increase A 2 by a factor 4, the first ratio of branching ratios increase 
to ~ 10 -3 without changing the other ratio, while still passing all the experimental cuts. 
In this way, it is clear that by measuring the branching ratios of the neutralinos we get 
information on the parameters of the model. 

The discussion above suggests that the observation of events coming from processes 
like pp -> x°iX°i -> WWer (at the LHC) or e+e" -> x°iXi -> WWer (at the NLC) would 
make possible to measure parameters relevant for neutrino physics. 

We use CompHEP 4.4 |24] to calculate the production cross sections a(pp — > X1X1) 
(LHC) and a(e+e~ -> X1X1) ( NL C at yfs = 500 GeV) at leading order. The relevant 
Feynman diagrams for the LHC are shown in Fig. |U1 For the SPS1 mSugra benchmark 
we obtain: 

<?{PP-+Xlx\) = 9.8xl0- 3 pb 
o-{e+e~ -> xSx?) = 0.27 pb (34) 
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The cross sections of the whole processes were calculated multiplying the production cross 
sections by the branching ratios B(xi ~^ W + e~) and B(xi ~^ W + t~). Their values, for 
the set of parameters we have chosen, are: 

B(x°i -> W + e-) = 1.9 x 10~ 2 

B(Xi -»• W + t-) = 5.9 x 10" 2 (35) 

The complete cross sections are: 

a(pp -> x°iX°i -> W" + W" + e-r-) = 1.1 x 10~ 5 pb 
a(e + e- -> -> W+W+e-T - ) = 3.0 x 10~ 4 pb (36) 

On the other hand, the main source of background comes from the production of four W's 
with two of them decaying leptonically. We calculated those processes using CompHEP 
and we found: 

a(pp -> WWWW -> iy + W^ + e-r-z7 e z7 T ) = 6.5 x 10~ 6 pb 
a(e + e- -> WWW -> W/ + W/ + e~r-z7 e z7 r ) = 1.6 x 10~ 6 pb (37) 

Assuming a luminosity of 10 5 pb/year at, both, the LHC and the NLC we expect 1 
signal event per year at the LHC and 30 signal events per year. Nevertheless we are not 
interested on the charge of the final leptons, we only require that one lepton belongs to 
the first family and the other to the third one, so the total number of signal events are 
obtained by multiplying the above results by four. 

We remark that while the background is small in both cases, the number of signal 
events at the LHC is also small and a more detailed analysis is required. On the other 
hand the NLC appears as very auspicious environment for studying this model. 



8 Conclusions 

We have re-examined the possibility of generating neutrino masses and mixing angles 
in Supergravity with bilinear R-Parity violation. We found solutions with a relatively 
large value of tan /3, such that one-loop contributions to the neutralino mass matrix are 
as important as tree-level contributions. The heaviest neutrino mass is still generated 
mainly at tree-level, but the other two masses and the three mixing angles are strongly 
affected by loops. In particular, the tree level approximations for the mixing angles give 
completely erroneous results. 

We concentrate our study on a texture for the neutrino mass matrix which is common 
among our solutions, and on one particular solution corresponding to this texture. The 
atmospheric mixing is nearly maximal, and the deviation of the parameter tan 2 #23 from 
unity is related to the smallness of the ratio between the solar and atmospheric mass 
scales Am 2 sol / Aml tm . In addition, the solar and reactor angles are both small because 
of the small parameter A, which in turn is small because Ax/A 3 < I. Nevertheless, the 
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reactor angle is much smaller than the solar angle because the second neutrino mass is 
much larger than the third one. 

We have shown how the neutrino observables depend on the BRpV parameters 6j and 
Aj, and this dependency can be understood in terms of simple approximations in terms 
of parameters A, B, and C, where all the complication of the one-loop contributions is 
concentrated. The dependency on q and Aj is strong, and it is not clear a priori that a 
solution is for granted, due to the increasing precision of the measurements of the neutrino 
observables. It is shown also how these observables depend on the Sugra parameters, 
namely the universal scalar mass mo and the universal gaugino mass Mxii- F° r the given 
by the values of 6j and Aj, solutions lie in a narrow strip in the plane m — M1/2, where 
the gaugino mass is strongly restricted by the solar and atmospheric mass scales, and the 
scalar mass by the atmospheric angle and mass scale. 

Finally, we showed how the decay rates of the neutralino depend directly on some of 
the parameters 6j and Aj. In fact T(xi ~^ We) and T(xi ~ > W/i) depend only on Ai 
and A 2 respectively, while T(xi ~^ Wt) depends on both A 3 and 63. Measurements on 
branching ratios of the LSP can therefore give important information on the parameters 
of the model. We estimated that a few events with e ± r ± in the final state can be observed 
at the LHC and about a hundred at the LC, indicating that a measurement of the decay 
rates is possible at the LC. A more detailed analysis is necessary to estimate the expected 
precision of these measurements. 
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